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We use a simple, instanton motivated, nonlocal chiral quark model to 
calculate pion Generalized Distribution Amplitudes (GDAs). The nonlo- 
cality appears due to the momentum dependence of the constituent quark 
mass, which we take in a form of a generalized dipole formula. With this 
choice all calculations can be performed directly in the Minkowski space 
and the sensitivity to the shape of the cutoff function can be studied. We 
demonstrate that the model fulfills soft pion theorems both for chirally even 
and chirally odd GDAs. The latter one cannot be derived by the methods 
of current algebra. Whenever possible we compare our results with the 
existing data. This can be done for the pion electromagnetic form factor 
and quark distributions as measured in the proton-pion scattering. 
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1. Introduction 

In this work we calculate two pion distribution amplitudes (27rDAs) 
and pion skewed parton distributions (vrSPDs) in the instanton motivated, 
effective chiral quark model. 27rDAs and vrSPDs are defined as Fourier 
transforms of matrix elements of certain light-cone operators, taken between 
the pion states. 27rDAs appear in the amplitude for the process 7*7 — > vrvr, 
if the virtuality of the photon is much larger then the squared invariant mass 
of the two pion system [I]. 27rDAs describe the transition from partons to 
two pions in a final state. They are also related to the pion electromagnetic 
form factor in the time like region and pion electromagnetic radius. 

The process 7*7 — > vrvr can be related by crossing symmetry to virtual 
Compton scattering j*7r — > -fTT. This process factorizes into a hard photon- 
parton scattering and vrSPDs In the forward limit the latter reduce 

to the quark densities which are measurable in the tt— p Drell-Yan process 
or prompt photon production Finally certain integral of vrSPDs gives 

pion electromagnetic form factor in the space like region. 

Pions are the simplest hadronic states being qq pairs and Goldstone 
bosons of spontaneously broken chiral symmetry at the same time. There- 
fore their properties may be calculated with little dynamical input, relying 
on their chiral structure and chiral symmetry breaking. In this work we 
use the instanton motivated effective chiral quark model [7] with nonlocal 
interactions. This model has been successfully applied to the calculation of 
the leading twist pion distribution amplitude 0-^0], 27rDAs ^J-^Sl and 
pion SPDs |12l I13j . The main ingredient of the model is the momentum 
dependence of the constituent quark mass which regularizes certain, other- 
wise divergent, integrals. To make the calculation feasible this momentum 
dependence has been taken in the form |9 : 

M is the constituent quark mass at zero momentum. Its value, obtained 
from the gap equation, is approximately 350 MeV. Quantities n and A = 
A(n; M) are model parameters. As explained below we expect the model 
to be roughly independent of n, if the value of A is properly adjusted. The 
formula maybe thought to be instanton motivated in a sense, that when 
continued to Euclidean momenta (A;^ —k\)-, it reproduces reasonably 
well [Sj the momentum dependence calculated explicitly |141 115j in the in- 
stanton model of the QCD vacuum: 

FinsiikE) = 2zMz)K,{z) - h{z)Ko{z)] - 2h{z)K^{z), (2) 

where z = pkE/2, with p = (600 MeV)""*^ being an average instanton size. 
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The calculations presented here may be viewed as an extension of the 
previous results of Refs.JJ-^31 and were partially reported in ^StiQ. Our 
motivations are both theoretical and phenomenological. First, the asymp- 
totics of Eq.Q suggests n = 3/2 in Eq.(^, however, taking the non-integer 
n results in additional difficulties. In the previous works [7| |S1 Illj-jl3j the 
value n = 1 was assumed. Therefore it should be checked that the results do 
not depend strongly on n. This is indeed true for the pion distribution am- 
plitude, as has been shown in _9". Secondly in [71IH1^2-J3| the momentum 
dependence of the constituent mass in the quark propagators was neglected. 
With the method developed in |9_ we can calculate generalized parton dis- 
tributions for arbitrary n, keeping track of momentum dependence of the 
constituent quark mass both in numerators and denominators even for the 
convergent quantities. This is in fact necessary if one insists on the soft pion 
theorems which would be otherwise not fulfilled. Apart from a well known 
soft pion theorem for a chirally even 27rDA , we derive (in the framework 
of the present model) a soft pion theorem for a chirally odd 27rDA. 

On the phenomenological side, apart from the calculation of the pion 
generalized distribution amplitudes (GDAs for short), i.e. 27rDAs and 
vrSPDs themselves, we calculate pion electromagnetic radius, electromag- 
netic form factor and quark distributions, and compare them with the ex- 
isting data. This comparison although not bad, is not entirely satisfactory. 
We shall discuss our results in Sect|71 We start with the short description of 
the model in Sect|2Then, in Sect|51 we proceed with definitions of 27rDAs 
and vrSPDs and their general properties. In Sect|l] we sketch the calcu- 
lations and present numerical results. In Sect Owe demonstrate soft pion 
theorems. Finally, in SectEl we discuss pion quark densities calculated in 
the model and compare them with data and other theoretical calculations. 
Technical details can be found in the Appendix. 



For two quark fiavors {u and d) the effective model we use to calculate 
pion GDAs is given by the action (in momentum space): 



2. Effective chiral quark model 



I 



(2^ 



I 



d'^p d'^k 
(27r)4 (2^ 



^ ^ (p) {p-k) \/Mfc V( fc) , 

(3) 



where 
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The matrix 

[/"^s (x) = exp 



— 7 T TT ( 



l + -^7Vvr«(x)-^^'^(xK(x) + 



(4) 

(in coordinate space) gives the interaction between quarks and pions. There 
is no kinetic term for the pions and the pion field is treated as an external 
source. = 93 MeV is the pion decay constant, are Pauli matrices 
and the pion T3 eigen states are: vr^ = vr'^, vr"*" = — (vr^ + ivr^) /\/2, vr" = 

(tt^ — iTT^) /\/2. The momentum dependence of the quark constituent mass 
is given by Eq.(^. The A parameter is fixed once for all by adjusting the 
value of Ft^{A) calculated in the effective model © to its physical value 

= 93 MeV. This has been done in jUj. For example for M = 350 MeV 
and n = 1 we have obtained A = 1157 MeV. 

Neither M nor A should be identified with the normalization scale Qq of 
the quantities calculated in the model. The precise definition of Qo is only 
possible within QCD and in all effective models one can use only qualitative 
order of magnitude arguments to estimate Qq. Arguments can be given that 
the characteristic scale of the instanton model is of the order of the inverse 
instanton size 1/p = 600 MeV, or so. We shall come back to this question 
in Sect El 

3. Definitions and general properties of GDAs 

We define two light-like vectors: = (1,0,0,-1) and n'^ = (1,0,0,1). 
These vectors define plus and minus components of any four-vector: = 
+ = n ■ k and k~ = k^ — k^ = n ■ k. 

3.1. Two pion distribution amplitudes 
We take the definitions of pion generalized distribution amplitudes from 

[HIS]. 

— exp{-iuXn-P)(^iT''{pi)TT\p2) i^^' iXn)T^^Jf (0) 0^ 
= 6'^' 5ff' ^i=\u, V, s) + ie'^'^ {r^yf ^=\u, v, s), (5) 

where P = pi + P2 is the total momentum of the two pion system. For the 
leading twist (chirally even) 27rDAs we have: 

r = ^, (6) 

whereas for the chirally odd 27rDAs 

r = in'^P'^a^f,. (7) 
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Note that with the definition Q, chirahy odd 27rDAs have the dimension 
of mass. 27rDAs depend on the fohowing kinematical variables: squared 
invariant mass of the two pions, s = P^; the longitudinal momentum frac- 
tion carried by quark with respect to the total longitudinal momentum, 
u = k~^/P~^ and the longitudinal momentum fraction carried by one of the 
pions with respect to the total longitudinal momentum, 
27rDAs have the following symmetries jl2j : 

^fiu, V, s) = -$^=°(1 - u, V, s) = $^=°(u, 1 - V, s), (8) 



^-'{u, V, s) = ^-'{1 - n, V, s) = -^-'{u, 1 - V, s). (9) 

The isovector chirally even 27rDA is normalized to the pion electromagnetic 
form factor in the time-like region: 

1 

I du^f{u,v,s) = {l-2v)Fns). (10) 



From its s dependence the pion electromagnetic radius can be evaluated 

(11) 



J pern / \ 



The normalization condition for the isoscalar 27rDA gives the fraction of a 
pion momentum carried by the quarks, M2 ■ 



du{2u - l)^{-^{u,v,s = 0) = -2v{l - v)M^''\ (12) 



It is convenient do expand the 27rDAs in the basis of the eigen functions 
of the ERBL equation QSJO^l (Gegenbauer polynomials Cn'^{2u — 1)) and 
partial waves of scattered pions (Legendre polynomials P;(l — 2v)): 

00 n+l 

<^Uu, V, s) = 6u{l -u)Y,Yl bUs)CII\2u - mil - 2v). (13) 

n=0 /=0 

Because of the symmetry properties ((HI El, the sum in Ea. (|13|l goes over 
odd (even) n and even (odd) / for isoscalar (isovector) 27rDA. 

The expansion coefficients Bj^^ are renormalized multiplicatively: 

bUs;Q) = bUs-M i^j^y , (14) 
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with 7„ being anomalous dimensions [SHI- Prom and (|14|) it is easy to 
read the asymptotic form of 27rDAs: 

^r^-'='iu,v,s) = (15) 

and 

<i>^'^-'=\u,v,s) = 6uil - u){l - 2v)bI,THs). (16) 

This asymtotic form is plotted in Fig^ For chirally even 27rDA, from H1U|) . 
we have B^^^s) = F^'^is). 

s = GeV2 




Fig. 1. Asymptotic form the isovector 27rDA as given by Eq.Cni. 

In Ref.jT2] the following soft pion theorem for the chirally even 27rDAs 
has been demonstrated. In the chiral limit (myr = 0), in the case of one of 
the pions momenta going to zero (w(l — ?;) ^ 0) we get 

^f{u, v = 0,s = 0) = ^f{u, v = l,s = 0)=0 (17) 

and 

^=\u, v = 0,s = 0) = -^=Hu, v = l,s = 0) = </>^^(n), (18) 

where (j)^^{u) is the axial- vector (leading twist) pion distribution amplitude: 

1 

(0|d(z)7^75n(-z)K+(p)) = iV2F^p^J dn e*(2.-i)..p^Ay ^^g^ 



We shall come back to this point in Sect|SJ 
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3.2. Pion skewed distributions 
For vrSPDs we have (for a review see Ref.{4,): 



— exp {iXXn ■ p) ( n^ip) 

ZTT \ 



f 



An 



An 



n'^ip) 



Sab^ff H^=^{X,i,t) + ie''^^ {t^Y^' H^=\X,^,t), 



(20) 



where p = {p + p')/2 is the average pion momentum. SPDs depend on the 
following kinematical variables: the asymmetry parameter ^ = — A'*"/(2p"'"), 
where A = p' — p is the four-momentum transfer; t = and the variable 
X defining the longitudinal momenta of the struck and scattered quarks in 
DVCS, (X + i)p+ and {X - ^)p+ respectively 

The pion SPDs obey the following symmetry properties jl3j: 



H'='>iX,C,t) = H'=^iX,-C,t) = -H'=^{-X,C,t) 



tI=0 



1=0, 



(21) 



H'=\X, t) = H'=\X, -C, t) = H'=\-X, e, t). (22) 

In the forward limit the SPDs reduce to the usual quark and antiquark 
distributions: 

H'=\X, e = 0, t = 0) = ^ [@{X)qs{X) - e{-X)q,{-X)] , (23) 

H'=\X, e = 0, t = 0) = i [e{X)q,{X) + e{-X)q,{-X)] , (24) 

where qs{X), qv{X) are singlet (quark plus antiquark) and valence (quark 
minus antiquark) distributions. For example for tt"*" we have: 

2H^-^iX,C = 0,t = 0) = \ (25) 

^ ' \ -(F {-X)-d'' {-X) for X<0 ^ ' 

and 

2/7-^(X,C = 0,t = 0) = | (26) 

^ ' ' \ -d^ (-X)+d^ (-X) for X<0. ^ ' 

In the present model we get u^^ {X) = d^^ {X) = 0. 

The normalization conditions for SPDs are analogous to the 27rDAs case 
(HnHnmil). in particular 

+1 

j dXH'=\X,^,t)=F^'^{t), (27) 

-1 
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Fig. 2. Diagrams contributing to the matrix element ^7r''(pi)7r''(p2) ^ip^ (zi)rtl;^ {z2)^ 
Eq.Q. Diagrams a) and b) contribute both to isoscalar and isovector 27rDAs; diagram 
c) contributes only to the isoscalar 27rDA. 



and 



// em\2\ r. ^-^n (^) 



dt 



+1 

J dXXH'=\X,i,t = 0) = ^(1 - e)M, 



2\ ;i/rW^ 



(28) 



(29) 



In analogy to the soft pion theorem for 27rDAs, in the case of ttSPD we have 



H^=^{X,i = l,t = 0) = 0. 



(30) 



4. Pion generalized distribution amplitudes in the effective model 

^.1. Analytical results 

For the matrix elements entering the definitions of 27rDAs and SPDs, 
Eqs. (|5l2()j) . in the effective model described in SectEl we obtain: 



^S^^ff [Ti{k - P,k)+T2{k- P,k- P2, k)+T2{k- P,k-pi, k)] 
^■^abc ^^cyf ^k-P,k- P2, k)-T2{k-P,k-pi, A;)]} , (31) 
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f,Z2 „ f',Zi 



f,Z2 ^ flZi 
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Fig. 3. Diagrams contributing to the matrix element ^7r''(p') |?/)^ (zi)r^''^(z2) | 7r"(p) 
Eq. 121)1 . Diagrams a) and b) contribute both to isoscalar and isovector SPDs; diagram 
c) contributes only to the isoscalar SPD. 



and 



A A 

Ti[k ,A;H 

^ * 2 2 



A _ A 
+T2 I k- —,k-p,k + — 



+ T2 k 



— ,k + p,k + — 



where 



A A\ / A A 

T2[k--,k-p,k + -\ - T^ik- -,k + p,k + - 



Tl{q,k)=Tv 



i-Mq + ie 



1 



%ir,q,k) = Tr 



-7 



5 



-7 



(32) 



(33) 
(34) 



f - Mr + ie' i- Mq + ie' ft- Mk + i€_ 

For the leading twist chiraUy even distributions F = fi. 

Equation ()32() is obtained, due to the crossing symmetry, from by 
the exchange pi — > —p and p2 — > p' (that is P — > A). The diagrams 
contributing to the matrix elements ()31() and ()32p are shown in figures [^J 
and 121 respectively. The diagrams a) and b) contribute both to the isoscalar 
and isovector GDAs, while the diagram c) contributes only to the isoscalar 
GDAs. 

In the case of 27rDAs we will work in the reference frame defined by 
Pi = 0. In this frame we find: 



P 



(35) 
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Pi = {vP+, (1 - v)^, p^) , P2 = ((1 - v)P+, v^, -p±) , (36) 

pi = ^(1 — v)s — m^. (37) 

From p^> it follows 

1 _ yr^^ <2v<l + ^1-^, s > Ami (38) 

In the case of pion SPDs we will work in the reference frame defined by 
p_i_ = 0. In this frame we find: 

Al = -t- fi-t + 4ml) = (1 - - ^rnl 

From > it follows 

In what follows we take the chiral limit, m-,^ = 0. In this limit, from (|38() 
and (gll), we get < -u < 1 and -1 < ^ < 1. 

Inserting Eqs.(|nil) and (EH into Eqs.® and we find A;+ = nP+ and 
fe+ = Xp'^ respectively. 

The method of evaluating dk~ integral, taking the full care of the mo- 
mentum mass dependence, has been given in 0. To evaluate dk~ integral 
we have to find the poles in the complex k~ plane. It is important to note 
that the poles come only from momentum dependence in the denominators 
of Eas. ()33l34() . Indeed, for each quark line with the momentum k, in the 
case Mk — > oo we have at most in the numerator (if the line is coupled 
to pion lines at both ends) and in the denominator. This means that 
the position of the poles is given by the zeros of denominator, that is by the 
solutions of the equation 



(39) 

(40) 
(41) 
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This equation is equivalent to 

4n+l 

G{z) = ^^"+1 + -r^ =J{{z-Zi), (44) 

i=l 

with z = k'^ / h? -l + ie and = M^/A^. For / (or finite A) equation 
(|44|) has 4n + 1 nondegenerate solutions. In general case 4n of them can be 
complex and the care must be taken about the integration contour in the 
complex k~ plane. Because of the imaginary part of the Zj's, the poles in 
the complex k~ plane can drift across Re/c" axis. In this case the contour 
has to be modified in such a way that the poles are not allowed to cross it. 
This follows from the analyticity of the integrals in the A parameter and 
ensures the vanishing of GDAs in the kinematically forbidden regions. We 
get nonvanishing results for 27rDAs and SPDs for < n < 1 and — 1 < X < 1 
respectively. 

s = GeV^ s = 1 GeV2 




Fig. 4. Chirally even isoscalar (upper row) and isovector (lower row) 27rDA's for M = 
350 MeV and n = 1 for two different invariant masses of outgoing pions' momenta. Note 
that for small s and v = Q the isovector 27rDA resembles - as function of u - the axial- 
vector one pion distribution amplitude in agreement with the soft pion theorem II18I63II 

After evaluating dk~ integral the dPk±_ integral has to be treated nu- 
merically. The integral over angular dependence in the transverse plane 
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can be in principle evaluated analytically (using residue technique), or in 
other words an exact algorithm for its evaluation can be given. In this way 
we are left with the numerical integration in only one variable, which is an 
easy task to do. The technical details of the calculations are presented in 
Appendix El 

s = GeV2 s = 1 GeV2 




Fig. 5. Chirally odd isoscalar (upper row) and isovector (lower row) 27rDA's for M = 
350 MeV and n — 1 for two different invariant masses of outgoing pions' momenta. Note 
that for small s and i; = the shape in u of the isoscalar 27rDA resembles the derivative 
(straight line) of the pseudo-tensor one pion distribution amplitude in agreement with the 
soft pion theorem l|68|l . 

4-2. Numerical results 

Our results for 27rDAs and SPDs are presented in Figs. |11|71 Although 
strictly speaking the present model is valid only for small pion momenta, s, 
we show also results for s as large as 1 GeV in order to trace the trends of 
the change of shape. 

For small s and one pion momentum equal zero {i.e. v = 0) the shape 
of the chirally even isovector 27rDA resembles - as function oi u - the axial- 
vector one pion distribution amplitude in agreement with the soft pion theo- 
rem (|18I63() . Similarly, chirally odd isoscalar 27rDA resembles the derivative 
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" ^ isoscalar, chirally odd 2jt distr. 
/ M=350MeV, 

' ''""■x^ J=0.25 GeV^ 

„=2 •'='>'25. 



«=1 



0.0 0.2 0.4 0.6 0.8 



0.5- 
0.4- 
0.3 
0.2 
0.1 



n=Hi\ 



isovector, chirally odd 
// 2jT distribution: 
/,' M=350 MeV, 
V .5=0.25 GeV\ 

v=0.25. 



0.0 0.2 



0.4 0.6 



Fig. 6. Isoscalar (left) and isovector (right) chirally odd 27rDA in the chiral limit for 
M = 350 MeV, n = 1 (dashed-dotted) and 2 (dashed line), s = 0.25 GeV^, at u = 0.25. 



(straight line) of the pseudo-tensor one pion distribution amphtude, also in 
agreement with the soft pion theorem (|68j) which we shall discuss in detail 
in Sectini The functions that we obtain, obey correct symmetry properties 
(ISl El) and . The asymptotic form for the isovector 27rDA is plotted 

in Fig^ and we see that it has a different shape from the model prediction 
of FigH 

As in the case of the pion distribution amplitude |3] the n dependence 
(see Eq.(^) of our results is weak. This is depicted in FigEl where we plot 
chirally odd 27rDA's for M = 350 MeV and n = 1 and 2, s = 0.25 GeV^, at 
V = 0.25. One can see that there is almost no n dependence except for the 
end point behavior of the isoscalar distribution which vanishes at u = 0, 1 
for all n > 1. This kind of behavior has been already discussed in Ref.|10| 
in the context of the pseudo-scalar one pion DA. The n dependence of the 
chirally even 27rDA is similarly weak. 

Skewed pion distributions for M = 350 MeV and n = 1 are plotted in 
FigElfor two values of t = and —1 GeV^. Note that we only plot them for 
< ^ < 1, the remainder in region of —1 < ^ < is symmetric according 
to Eas. (|22l21j) . In the present model 

rri=^(X f-n f-n)-l ^'=°(^,C = 0,t = 0) for X>0 
H (X,^-0,t-0)-| _^,=0(^^^^0^^^0) fo^ x<Q 

and from (^11^ it follows that aJX) = qv{X). The section of i?^=°(X, i) 
as function of X is plotted in Fig|Hl In Fig|Hla we plot H^=°(X, ^ = 0.5, t = 
—0.25 GeV^) for M = 350 MeV and n = 1. The solid curve corresponds 
to the full result, whereas the dashed-dotted curve corresponds to diagrams 
a)-|-b) of FigEl and dashed line to diagram c). In the forward limit (t = 
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and ^ = 0, FigEJb) H^^^{X, ^ = 0, t = 0) corresponds to the quark densities 
in the pion (^5]) . We see httle n dependence. For comparison we also plot 
quark distributions obtained in the model with sharp cutoff in the transverse 
momentum plane which coincides with the result obtained in Refs.|2ll I22j. 
These quark distributions are understood to be at low normalization scale 
Qo and have to be evolved to some physical scale Q at which they are 
experimentally accessible. This will be done in SectEl 

Our skewed distributions do not exhibit factorization [21] and for ^ = 
they are quite similar in shape to the distributions obtained from local 
duality 



t = GeV^ t = -1 GeV^ 




Fig. 7. Isovector (upper row) and isoscalar (lower row) pion SPD, H'=^'''{X,^,t), in the 
chiral limit (m^ = 0) for M = 350 MeV; for i = and -1 GeV^ 

The values of F^'^(s) and Ff^^(t) evaluated by means of Fas. (1101271) do 
not depend (as they should) on v and respectively. They are depicted in 
FigEl together with the experimental data of Refs. |23]-[17j- We see that 
the space like pion form factor overshoots the experimental points. This 
discrepancy may be eventually cured by negative contribution from the pion 
loops. In the time like region we see that very soon the p resonance tail 
switches on. Of course vector mesons are not accounted for in the present 
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model. The similar buildup of the p tail can be seen in the vr — vr scattering 
data |2H|- One should note that strictly speaking our model can be used only 
for very low momentum transfers. For higher momenta QCD perturbative 
techniques should be used [231 

The values of the pion electromagnetic radius, obtained from Eas. ()ll() 
and H28() are consistent with each other, although too small. Numerically 
we get ((r^"*)^) = (0.54 fm)^, which is comparable with the value 

i^rT?) = 4^ = (0-58 fm)^ 

obtained within the instanton models. The experimental value |25j is, how- 
ever, bigger: {{r^)'^) = (0.66 fm)-^. This can be, however, cured by the 
pion loops, which are neglected in the present approach. 

5. Soft pion theorems 

Soft pion theorems relate 27rDAs for one pion momentum going to 
with the one pion light cone distribution amplitude. 

We recall the definitions of the axial-vector and pseudo-tensor one pion 
DAs: 

+ 00 

^^y{u) = -^ [ -e-''^'^-''>'^--^U0\d{Xn)fij,n{-Xn)\7r+{P)) (46) 
zV2F^ J TT 

— oo 

+ 00 



du ^ ' iV2{qq) J vr 

— oo 

X (0| d{Xn) in^^P^a^p-i^ n(-An) |7r+(P)) . (47) 

The matrix element entering the definitions 1)461 147(1 . calculated in the 
present model, gives: 

(0|d>i)7n(z2)|vr+(P)>=-^ j ^e^^^"^)--^- ^ A: - P), 

(48) 




Mk 5 \/^ 



- Mfc + iP iit-p)- Mk^p + iej 

^MkMk.p Tr (77' [-^ + MkW -n + Mk-P]) 
[k^ - Ml + ie] [{k - Py - Ml_p + ie] 



k-P 



(49) 
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Fig. 8. Isoscalar pion SPD, H'-°{X,^,t), in the chiral limit {m„ = 0) for M = 350 MeV; 
for t — —0.25 and GeV^ plotted as function of X for fixed ^. In Fig. a) dashed-dotted 
curve corresponds to diagrams a)+b) of Fig. 2, and dashed line to diagram c), whereas the 
solid curve corresponds to the sum. In the forward case (Fig.b)) we show results for n — 1 
(solid) and 5 (dashed-dotted line), as well as the distribution corresponding to constant 
M{k) Note that in Fig.b.) 2H'^° for X < is equal to -d{-X) distribution in tt+ 
while for X > it corresponds to u{X) . 



17 



Further calculations require an explicit form of 7. Substituting (|48j) into 
(gni) and (gZI) we get 




-t [GeV'] 



10- 



M=350MeV,n = l i § / 

i ' 



! / 



0.0 0.2 0.4 0.6 

."^ GeV 



0.8 1.0 



Fig. 9. Pion form factors in space like (upper fig.) and time like (lower fig.) regions. 
Theoretical curves correspond to M = 350 MeV and n = 1. Experimental data are taken 
from: dots triangles squares |27|. 
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and 

respectively. Similarly, in the case of 27rDAs, from © and 1)31^ we get 

[T^{k -P,k)+T2ik-P,k- P2, k)+T2ik-P,k-pi, k)] 

(52) 



and 



H=\u,v,s) =^ / -^6{n.{k-uP)) 



J {2^y 

[T2 {k-P,k- p2, k)-T2{k- P,k-pi, k)] , (53) 

with Ti and T2 defined in (|33|) and (|34|) . Soft pion theorems relate matrix 
elements (|48j) and (|3ip for one pion momentum, say p2 ^ and for different 
7 and r. For 7i we have: 

^n. o .^ _ V^fc-P^fc Tr (F [(^ - /^) + Mfc_p] + M^]) 

' ^" [{k-Pr-Mlp + ze\[k^-Mi + ie\ ' ^'^^ 

whereas for the two 72 's in the limit p2 ^ (that is v ^ 1) we get: 



TM^I^M, Tr (r [(^ M,_p]) 

' ' T(^^7^F^ic7^^^n^^^w^' 

r2(fe -P,k-P,k) = - ^^^_ _ ^^^^ _ ^^^^ . (56) 

5.7. Chirally even 2ttDA and axial-vector wave function 

Consider first the soft pion theorem discussed already at the end of 
Sect 13. 11 To this end let us take 

7 = ?^75- (57) 

Then 

r(kh-p\\ - ^^k^k-p [-^k-p TV {m + TV (^(^ - m 

^' ^^lr=^5 [k2_Ml + ie][{k-Pf-Ml_p + ie] 

(58) 
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For two pions we have to take 

T = fi. (59) 

Then in the soft Hmit pi = P, p2 = (this means v = 1) the isovector 
combination in Eg . (1531) reads 

T2{k - P, k, k\^^ - T2{k -P,k-P, k)\^^^ 

[{k - Py - Ml_p + ie] [k^ - Ml + ie] 
= -T{k-PM,=^,, (60) 
Similarly, for the sum entering the isoscalar combination in Ea. ()31() we have: 
Ti{k - P, k)\^^^ + T2ik - P, k, k)\^^^ + T2{k -P,k-P, k)\^^^ = 0. (61) 

Hence the soft pion theorem for the chirally even 27rDA takes the following 
form 

H^"^eveniu,V = l,S = 0)=0, (62) 
^itxevenin, V = I, S = 0) = -0^^^(n). (63) 

5.2. Chirally odd 271 DA and pseudo-tensor wave function 

Let us now consider chirally odd 27rDA in the soft pion limit. To this 
end let us take 

^ = m°PV/375 = ^[/',?i]75 (64) 

whose matrix element defines the derivative of pseudo-tensor one pion dis- 
tribution amplitude. Then 

T(i.ie-PM 1 ^MkMk-pT,[mm-r) 

M/^,^ ^J^=|[^,^]^3- 2[k^-Ml + ie][{k-Pf-Ml_p + teY^ ' 



For 2 pions let us take 



Here only Ti survives 



r = m°PV/3 = ^[f,?i]. (66) 



^il^; ^, ''^lr=i[AT^] 2 [(A: - Pf - Ml_p + ie] [k^ - Ml + te] 



nk,k-P)\^^,j^.^ (67) 
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while T2 = 0. Therefore the soft pion theorem relates the isoscalar part of 
the chirally odd 27rDA to the derivative of pseudo-tensor one pion amplitude: 

H='^odAu,v = l,s = 0) = J^±<l,PJ{u), (68) 

^itxoM{n,v = l,s = Q) = Q (69) 

This result cannot be derived by current algebra. Since the pseudo-tensor 
one pion distribution amplitude is close to the asymptotic form 6m (1 — u) 
[TU] therefore we expect ^^^^^odA'^-'''^ = 0,s = 0) ~ 1 — 2u. That this is 
indeed the case can be seen from FigEJa. We recall that with our definitions 
© and (O the chirally odd 27rDAs have the dimension of mass. 



6. Pion structure function 

As explained in Sect 13.11 skewed pion distributions are in the forward 
limit directly related to the parton distributions inside the pion. For exam- 
ple, as explained after Eas. (|25l26|) and as seen from Fig|Hlb 

n^+(x) = dr^{x) = 2H^=\x,^ = 0,t = 0). (70) 

On the other hand one can relate pion quark distributions to the wave 
functions for all Fock states and polarizations [30]. For 7r~^ we have 



Fock states ^ 



2(27r)' 

V^^^(x,A:t)|' + ... (71) 



where we have saturated the sum over the different Fock components and 
polarizations by the axial-vector wave function. The function ij)^^ is (for 
A'^c = 3) normalized in the following way j3Uj : 



/ 



1 

'2 





It is clear that the normalization condition for m" 

1 

cixu^'(x) = l (73) 



1 
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and normalization condition (|72j) are in general incompatible since other 
components denoted by dots in Ea.H71|) are of importance. For example in 
a model with a constant M one obtains for ifj'^'^ 7\: 



'AV\ 



e(n(l-n)) 



2^/3 M2 



(74) 



and the normalization condition (|72|) is achieved by imposing a cutoff A on 
the transverse momentum integration dk'^: 



m ■ 



(2vr)^ 

It is now straightforward to calculate u'" 



^+(x) =e(a;(l-a;)) 



6 

(^^ 



M2 



1 — exp 



(75) 



(76) 



The normalization constant is, as expected, smaller than 1, however, for 
M = 350 MeV we get 0.9 - a fairly satisfactory result for such a simplistic 
model. Of course the properly defined quark distribution is also properly 
normalized i21, 22_. 

In the present model |H] 



2V3M^ 



+ 1 + ZiX + Zk{l - X 



(77) 



and the normalization condition for M = 350 MeV and n = 2 gives 0.88 
instead of 1. The shape is also different from the properly normalized result 
obtained with the help of Eq. (|70j) . however, it can be explicitly seen that 
for large x both definitions converge, as they should [HJ. This is depicted 

in Figirni 

One of the major problems of the effective models like the one considered 
here, is the normalization scale Qq at which the model is defined. This is 
crucial shortcoming as far as the comparison with the experimental data is 
concerned. It is argued that the relevant scale for the instanton motivated 
models is of the order of the inverse instanton size 1/p i.e. approximately 
600 MeV. The precise definition of Qq is only possible within QCD and in all 
effective models one can use only qualitative order of magnitude arguments 
to estimate Qq. A more practical way to determine Qq was discussed in 
Ref. jlO| where we associated Qo with the transverse integration cutoff Kt 
which was of the order of 760 < Kt < 1100 MeV. 
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M=350 MeV, n=2 




Fig. 10. Comparison of the u quark distribution in vr"'" defined by Eg. 1)70^ (solid line) and 
the "naive" one defined by Eq.||7IJ (dashed line) for M = 350 MeV and n = 2. 

On the other hand in Refs.[22I it was argued that the pertinent model 
scale may be as small as 350 MeV. This estimate is based on the requirement 
that the total momentum carried by the quarks equals to the one measured 
experimentally at Q = 4 GeV^, i.e. 47% . In that way the initial evolution 
scale Qo can be adjusted. It is, however, problematic whether one can use 
QCD evolution equations at such low Qq. 

Since, as will be discussed in SectQ, our model has problems with the 
momentum sum rule, therefore we cannot use the above prescription to fix 
Qq. In Figs llll we simply show the shape of the valence, sea and gluon 
distributions calculated in the model and evolved (in the leading log ap- 
proximation) to the scale = 4 GeV'^ assuming Qq = 450 and 350 MeV. 
As initial conditions we take valence quark distributions as calculated in the 
model with sea and gluon distribution equal to zero at the initial scale Qo- 
Therefore both sea quarks and gluons are generated dynamically during the 
evolution. In Figs llll we also show "experimental data" as extracted from 
the pion-proton Drell-Yan and direct photon production [SJ [S] . For compar- 
ison we also show results of Ref.j21| with a constant quark distribution at 
initial scale. 

It can be seen from Figs llll that our quark distribution differs from the 
distributions extracted from the data. However, two existing experimen- 
tal parameterizations are not compatible. Interestingly, the constant initial 
quark distribution after evolving to Q'^ = 4 GeV^ fits very well param- 
eterization of Ref.[5]. Unfortunately for the sea, both constant and our 
distributions do not follow the experimental parameterization. This sug- 
gests that rather than compare model results for pion parton distributions 
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with the ones extracted from the data it is perhaps more appropriate to 
calculate the cross section itself and compare directly with the data. 

There have been several other calculations of the pion structure func- 
tion in the literature. A result similar to ours was obtained in Refs.|S2] 
where the authors calculate explicitly a hand-bag diagram in a local NJL 
(nonbosonized) model in the Bjorken limit. They, however, introduce an 
X dependent cutoff in order to regularize the divergent integrals and get 
proper behavior of the valence quark distribution in the large x limit. Their 
result is in agreement with a similar calculation of Ref . j^Sl • In an approach 
based on Ward-Takahashi identities, which is in fact equivalent to the local 
NJL model with a sharp momentum cutoff, the valence quark distribution is 
equal 1 over the whole range of x j22] . This very simple quark distribution 
is properly normalized, in a sense that the quark number is 1 and its total 
momentum is 1/2. The vanishing of q{x) for x ^ 1 is achieved by DGLAP 
evolution. We have plotted the result of this evolution in Figllll 

Direct calculations in the instanton model |34j show phenomenologically 
quite similar behavior as ours. An advantage of Refs.pi] is that they use 
well defined currents with nonlocal pieces [SHI, whereas we use naive quark 
billinears. This is reflected in wrong normalization of the first moment of 
the quark distributions which at low scale should be 1 (for / dxx{u + d) ) as 
opposed to 0.93 what we get. An even larger mismatch has been reported 
in a similar model of Ref. j36j , where the remainder of the momentum was 
attributed to gluons and sea, which are, however, absent in our approach. 

Our quark distributions show little n dependence. This is depicted in 

Fig. na 

7. Discussion 

In this paper we have calculated the Generalized Distribution Ampli- 
tudes using the nonlocal, instanton motivated, chiral quark model. The 
nonlocality has been taken in the form of Eq.(^ which allowed us to perform 
all calculations in the Minkowski space. Introducing light cone integration 
variables allowed us to perform most of the calculations analytically. As a 
result we were left with the numerical integration in only one variable, dk^. 

We define the GDA's through the matrix elements of the nonlocal quark 
billinears (|5I2U|) . Although in QCD this is certainly a correct definition, one 
might envisage alternative definitions which would be more appropriate for 
the nonlocal effective models |211 134j . The reason is that in the limit when 
the quark operators are taken in the same point, quark billinears 1)5120^ do 
not correspond to the properly normalized Noether currents. This is be- 
cause in the nonlocal models Noether currents get additional pieces which 
restore Ward-Takahashi identities It is not clear how to generalize 
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Fig. 11. TT^ valence, sea and gluon momentum distributions. Upper panel: input dis- 
tribution as calculated in the model for M — 350 MeV and n = 1 (upper red dashed- 
dotted curve) and after evolution to = 4 GeV^ from Qo = 450 MeV (middle blue 
dashed-dotted curve) and 350 MeV (lower purple dashed-dotted curve). Green dashed 
line represents constant initial quark distribution 1221 . Experimental parameterizations 
are depicted by a yellow band [5] and orange line Lower panels: the same for sea 
quarks and gluons with input distributions equal zero. 



quark billinears ()5I2U|) . since such generalization is, in principle, process 
dependent. Therefore an alternative way would be to calculate the whole 
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Fig. 12. The n dependence of valence quark distributions for AI — 350 MeV and n = 1 
(solid) and 5 (dashed-dotted). Input density corresponds to the model result. Two other 
sets of curves marked by Qo — 450 and Qo = 350 MeV correspond to the distributions 
evolved to = 4 GeV^ assuming that the model scale was Qo- 

scattering amplitude directly in the effective model and then extract the 
quantities one is interested in by imposing certain kinematical constraints, 
like Bjorken limit for example. Although this procedure seems at the first 
sight attractive, there is a problem because the Bjorken limit requires large 
momentum transfer, whereas the effective models are defined at low mo- 
menta. 

In fact the precise determination of the normalization scale Qq at which 
the model is defined poses a serious problem. In the instanton model of 
the QCD vacuum that the pertinent energy scale is of the order of the 
inverse instanton size Qo = l/z^ ~ 600 MeV. However, it has been argued 
in Ref . j22j that Qq may be as low as 350 MeV. This estimate is based on 
the requirement that the valence quark distribution calculated in the model 
of Ref.j^ and evolved from Qo to Q = 2 GeV carry observed fraction of 
total momentum. Unfortunately we cannot apply the same procedure in our 
case, since the momentum sum rule is violated in our model. Prom both 
equations, (|12|) and (|29j) we find the momentum fraction carried by quarks 
to be 93%, independently of v and ^. This causes the problem, as in the 
model we use, the pion is built from constituent quarks (there are no gluons) 
so there is 7% of the pion momentum missing. A natural explanation of this 
mismatch is that we missed some terms which, in the limit where the quark 
operators are in the same point, would reduce quark billinears H5I2U() to the 
full nonlocal Noether currents |351 13 7j . 
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Appendix A 

Technical details of the calculation of the pion GDAs 

Inserting (|ST|) and into ((SJ and (^0)1 respectively we get 

^=^{u,v,s) = Ji{u,v,t) +J2{u,v,t) +J3{u,v,t), 
^i=\u,v,s) = J2iu,v,t) - Ji{u,v,t) (A.l) 

and 

H'=^{x, e, t) = ii{x, e, t) + i2{x, t) + is{x, t), 

H'=\X,^,t) =Ii{X,^,t) -l2{X,C,t). (A.2) 
Ji, J2, J3, ^1, ^2, stand for the integrals: 



uP+ 
(A.3) 



uP+ 
(A.4) 

(A.5) 



Ji{u,v,s) 



iNc 
2(2vr)4F2 



+00 

j d^k± j dk^T2{k- P,k-pi,k) 



J2{U,V,S) 



iNc 
2(2vr)4F2 



j d^k^ j dk^T2ik- P,k-p2,k) 



J3{U,V,S) 



iNc 
2(2^)4F2 



+00 

j d^k± j dk^Ti{k-P,k) 



k+ = uP^ 
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Ti(X,C,i) 



iNc 



4(27r)^F2 



cfk± / dk^ T2 k 



iNc 



(fki 



iNc 



4(27rrF2 



i-cx 



A _ A 

— ,k — p,k -\ 

2 ' 2 



A A 

dA:_T2 ( /^;--,A;+p,fc + - 



A A 

dk^Ti[k - —,k+ — 



k+ = Xp^ 
(A.6) 



k+ = 
(A.7) 



(Ai 



with 7i and T2 defined in ((331 OH- It is convenient to introduce the scaled 
variables: 



k^' 



A' A2 
and the notation: 

We introduce the factors 



P± 



t 



A ' 4A2 ^ 

u = 1 — U, V = 1 — V. 



4 = ^ (A.9) 



(A.IO) 



4rt+l 



k=l 



(A.ll) 



which we will use below. Zi are 4n + 1 roots of the Eq. (|44jl . The factors fi 
have a property: 



4n+l , 



for m < 4n, 
for ni = An. 



(A.12) 



This property is crucial for the convergence of the integrals (lA.SIIA.S)) . in 
analogy to the case of the pion distribution amplitude This property is 
true for any set of different 4n + 1 numbers, irrespectively to the fact that 
they are solutions of certain polynomial equation. 
It can be shown that 



J3(u, V, s) = Jsiu, s) = -Jsiu, s), 



(A.13) 
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Ji{u,v,s) = -J2{u,v,s), (A. 14) 



Ji{u,v,s) = J2{u,v,s) (A.15) 



and 



T2{-X,i,t) = -IiiX,C,t), M-X,i,t) = -l3iX,C,t), (A.16) 



l2iX,-^,t)=Ii{X,^,t), l3{X,-C,t)=l3iX,C,t). (A.17) 

The properties ()A.1,3IA.17|) . together with Eqs. ()A.1IA.2|) . provide the cor- 
rect symmetry properties for 27rDAs (0 and SPDs (EH 1121) 
For J^{u,s) we get analytical results: 



TV 4n+14n+l 

= 7:7^2^ "^^ "^fc (^^*^" ~ ^^fc") (-^^'^^ + l + UZi+ UZk) 

(A.18) 

If < ti < ?; then i72('U, v, s) reads 



r27r)3F2 l^^iJi ^i'^^i M 

fc=l 

C(o,i) n {A{u)\^ + B,k{u,v)\ + A{u)) 

k=l 

--■.T{u,v,s), (A.19) 
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where 



giu,v) = 



u 



2n 



-U I ^ I + UZi 



[U — V 



A{u)X'^ + bi{u, v)X + A{u) 
An 



2n 



2n 



unu)"^ + k\ + \ + uzi\ / A{u)\^ + bi{u,v)X + A{u) 



An 



-vuuuj'^ + {u + u — v) k\ + u{\ + Zi) + {u — u) \J k'^t]_ 



A2 + 1 



2A 



(A.20) 



If n < n < 1, then 



J2{U,V,S) = -J^{U,V,S), 



(A.21) 



with T{u,v,s) defined in HA.19() . The symbols A{u), bi{u,v), Bik{u,v) in 
Eq. (TOOl rOnl) stand for 



A{u) = u^J k\t 



2 ^2 



&j(n, v) = nr^ + nt'(n — v)u:'^ + v{k\ + 1) — (n — v)zi, 
Bik{u, v) = ut]_ + uv{u — v)uj^ + v{k\ + 1) — (n — v)zi + uz^. 

Because H^^^ and H^^^ are symmetric in ^ we can assume ^ > 0. Then 
is nonzero only if — ^ < < 1- Similarly X3(X, is nonzero 
only if < X < ^. 

For -(,<X<^, T^{X,i,t) and Ji(X,^,t) read 



C(0,1) 

jX + [(X + g) A z,]^" + (X - e) (^A^ + fe.A + Ay 

4n+l 

n [AA2 + i?ifcA + A] 

k=l 



(A.22) 
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■AT A/f2 °r 

x,(x,e,t) = (-ir^g5^/d(4) E/' 



c(o,i) n + B,feA + ^) (CA2 + AfeA + C)] 

k=l 



where 



(A.23) 

For ^ < X < 1, Xi{X,i,t) reads 



n 2 — 1 



Xi(X,C,i) = (- 



[(X - [C\^ + /,A + Cy (CA^ + ff,A + Cy ^(,) 

4n+l 1 ^ 'S) /! 

c(o,i) n [(CA2 + FifeA + C) (CA2 + GifcA + C)] 

fe=i 

(A.24) 

where 



(x-i)A ; 



2n 



+ ' (x-i)A ; (x-i)A ; "^'^ 



. (A.25) 



31 



The symbols A, bi, Bik, C, di, Dik, fi, Fik, gi, Gik in Eqs. ()A.22I - 1X25]) 



bi = 2( [{X' -f)r + Kl + 6l + l]+{^- X)zi, 
B,k = 2i [{X^ -f)r + Kl + 5l + l]+{i- X)zi + + X)zk, 

C={l-X)^^^, 

di = {X- 1) [(e + X){1 - Or + 6l] -{l + O (4 + 1) + - l)z^, 
Dik = (X-l) [(e + X)(l - Or + 6l] -(1+0 (4 + l)+{X-l)z^-{^+X)zk, 

fi = {X- 1) [iC + X){1- Or + dl] -{1+0 (4 + 1) - + X)zi, 
F,k = {X-l) [(e + X){1- Or + 4] -(1+0 (4 + l)-{^+X)zMX-l)zk, 

gi = {x- 1) [(e -x){i+ Or - 4] + (1 - e) (4 + 1) - (e - 

Gik = {X-l) [{i - X){1 + Or - 4] +(1-6 (4 + l)-{i-X)z,-{X-l)zk. 



stand for 
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